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We present a semiclassical explanation of the so-called Bohigas-Giannoni-Schmit conjecture which 
asserts universality of spectral fluctuations in chaotic dynamics. We work with a generating function 
whose semiclassical limit is determined by quadruplets of sets of periodic orbits. The asymptotic 
expansions of both the non-oscillatory and the oscillatory part of the universal spectral correlator are 
obtained. Borel summation of the series reproduces the exact correlator of random-matrix theory. 



PACS numbers: 05.45.Mt, 03.65.Sq 



Quantum spectra of individual chaotic systems can be 
phenomenologically described in terms of random-matrix 
theory (RMT) 0, . This universality - asserted by the 
celebrated Bohigas-Giannoni-Schmit conjecture (BGS) 
3] - is an empirical fact, supported by a huge body of 
experimental and numerical data. Proving its concep- 
tual origin remains one of the fundamental challenges in 
quantum or wave chaos. 

Spectral fluctuations are conveniently characterized 
in terms of the two-point correlation function, i?(e) = 
A 2 (p(E + ^)p{E - §§)) - 1, where p(E) is the energy- 
dependent density of states, A = l/(p) the mean level 
spacing and (■) denotes averaging over the energy E. Pre- 
dictions made by RMT for the two-point correlation func- 
tion are fully universal in that they depend only on the 
parameter e, and the fundamental symmetries of the sys- 
tem under consideration. Specifically, the complex repre- 
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sentation R(e) = lim 7 _o R- e C( e+ ) where e d 
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and 



C*(e+) = £ ?r (ttG + (E + eA/2w)ix G~(E 
is employed, with G ± (E) = (E ± ijA/2ir - H)' 1 
H the Hamiltonian. The Wigner-Dyson unitary and or- 
thogonal symmetry classes (the only ones to be consid- 
ered here) of RMT afford the asymptotic series 
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In either case, C(e + 
(power series in \ je 



is a sum of a non-oscillatory part 
and an oscillatory one (e 2le+ times 
a series m 1/ e + ). Rcsumming the series by Borel summa- 
tion techniques and extrapolating to small positive values 
of e, one obtains R(e) + 1 oc e 13 for e — > 0, a signature of 
the level repulsion symptomatic for chaos ((3 — 1,2 for 
the case of orthogonal, unitary symmetry, resp.) 

The question to be addressed below is how to obtain 
the RMT prediction for a concrete chaotic (globally 
hyperbolic) quantum system. A step in this direction was 
recently made on the basis of Gutzwiller's semiclas- 
sical periodic-orbit theory p|. Gutzwiller represents the 
level density p(E) as a sum over periodic orbits, where- 
upon the function i?(e) becomes a sum over orbit pairs 



0, 0) Q • Relevant contributions to that double sum were 
shown Q to originate from orbit pairs which are identi- 
cal, mutually time-reversed, or differ only by connections 
in certain close self-encounters. By summing over all dis- 
tinct families of orbit pairs, the Fourier transform of R(e), 
the spectral form factor K(t), was found to coincide with 
the RMT prediction for times t = tTh smaller than the 
Heisenberg time Tjj = 27rS/A, the time needed to resolve 
the mean level spacing. The behavior of K(t) for r > 1, 
also known from RMT, was left unexplained. 

We now want to fill the gap left. As results we will 
obtain the full expression for the correlation function in 
the case of unitary symmetry, and an asymptotic 1/e- 
expansion amenable to Borel summation in the orthogo- 
nal case. The oscillatory term, Fourier transformed, then 
complements K(t) to its full form at r > 1. In many re- 
spects, our reasoning is inspired by the field theoretical 
formulation of RMT correlation functions , notably the 
existence of "anomalous saddle points" in the nonlinear 
cr-model ^(j . It also affords a new interpretation of ideas 
underlying the "bootstrapping" |Tlj . 

The basic idea of our approach is to consider represen- 
tations of C(e + ) different from the standard one in terms 
of the product of a single retarded and advanced Green 
function. We start from the generating function 
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where E^ b c d are energies in the vicinity of E defined 



= E + B c d A/2tt. From Z, the complex 
correlator can be accessed as 
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The two derivatives produce the product (trG + )(trG~) 
under the energy average. If we subsequently identify the 
energies "columnwise" (||): e\ = = e+, = = 
— e + , or "crosswise" ( x ) : e\ = e + , = — e+, = — e~, 
f^j = e~ , 7 — » the ratio of determinants approaches 
unity. The first representation for C(e + ) does not even 
require the limit 7 — > 0; it is widely used in RMT. The 
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second representation is crucial for us. Importantly, the 
semiclas steal approximation of either of these two exact 
representations misses contributions to Z, and therefore 
to C(e + ): the first representation yields only the non- 
oscillatory contributions, and the second (without "—5") 
only the oscillatory ones; adding both we will recover the 
universal two-point correlator. 

To see the emergence of these structures, let us repre- 
sent the determinants in as 

r r E i 1 

det(E% -H)~ x =exp| - / dEtrG + (E)j 
~ const x exp (inE^/A + ^ F a e^ s "^A , (4) 

a 

where the last line invokes the semiclassical expansion 
of the integrated Green function into a smooth (Weyl) 
average and a fluctuating (Gutzwiller) part; the latter is 
a sum of periodic orbits a with action S a and stability 
amplitude F a ; for simplicity, we assume the average level 
density 1/A to be constant; the periodic-orbit sum con- 
verges for ln\E\ large enough; the "const" in (0J comes 
from the lower limit of the energy integral and cancels 
from the ratio of determinants in Z . 

Expanding the exponential in (0J we get a sum over 
non-ordered sets of orbits A ("pseudo-orbits"), 
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A pseudo-orbit A may involve any number ua of compo- 
nent orbits (riA = pertains to the empty set which con- 
tributes unity to the sum); Fa is the product of the stabil- 
ity amplitudes and Sa the cumulative action of all com- 
ponent orbits. Expressing all four determinants in (0) 
similarly to llll- r >[) (e. g., using det(-E^ - H) = (det(E^ - 
H))*) and writing S(E + eA/27r) - S(E) + T(E)eA/2n 
(T is the period of an orbit, or the sum of periods in a 
pseudo-orbit) we approximate the generating function as 

Z ~ e 5 (e A- e s- e c+ e i)) (6) 

x (exp ( J^a F a e^ Sa(E)+l ^ e+A + £ b F 6 *e~* S * (K)_i *£ e5 
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Here, the mean density produces a phase factor 
ei( e A _e B _e c+ e r>). When representing the correlator as 
in 0, that phase factor turns into 1 and e 2le for the 
columnwise and crosswise identifications of energies, re- 
spectively. Indeed, then, we can recover cither the non- 
oscillatory or the oscillatory contributions to C(e + ). 

Another phase factor involves the difference AS = 
S A (E) - S B (E) + Sc(E) - S D {E) between the cumu- 
lative actions of (A, C) and (B,D). Due to this factor, 



systematic contributions in the limit h — > can arise only 
for quadruplets of pseudo-orbits whose action difference 
is of the order of h or smaller. 

The most basic of quadruplets have each of the com- 
ponent orbits of A and C repeated in cither B or D, 
such that AS = 0. These "diagonal quadruplets" may 
be summed by a lowest-order cumulant expansion: de- 
noting the periodic-orbit sums in the exponent of JfJJl by 
X, we may write (e x ) diag = exp{(A 2 ) diag /2}, wherein 

(-X" 2 ) diag contains only pairs of identical orbits. We ob- 
tain 



Z diag ~ exp(A 2 ) diag /2 
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Relying on ergodicity, the resulting sums over orbits may 
be evaluated by the sum rule of Hannay and Ozoria de 
Almeida [13, EJ^| 2 (-) » the lower limit 

of the integration is a minimal period To starting from 
which periodic orbits behave approximately ergodically. 
By this rule, e.g., the first sum turns into — ln(i(ej[ — 
+ cons t + 0(H). All four sums yield 
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with k = 1 for the unitary class. For the orthogonal class 
we must also consider pairs of mutually time reversed 
orbits. Therefore, each sum in JSJl must be multiplied by 
2 whereupon in the final result I© we have k = 2. 

Taking derivatives and columnwise identified energies, 
we recover the leading non-oscillatory contribution to the 
two-point correlator, (ie + ) _2 //3. Crosswise identified en- 
ergies yield the oscillatory contribution — c 2le (ie)~ 2 /2 for 
(3 = 2 (thus completely reproducing QJ), while for (3 = 1 
we get zero, i. e., no oscillatory term arises up to 0(e~ 2 ), 
again as in RMT. 

Going beyond the above level of approximation, we 
note that small phases may also arise from component 
orbits B and D differing from A and C in topology, but 
only weakly in action. By way of example, consider 
a pseudo-orbit A with just one component, while C is 
the empty set. Assume that A contains a so-called "2- 
encounter" , where two stretches of the orbit come close in 
phase space. This situation is depicted in Fig.QJi (where 
the phase-space encounter contains a crossing in config- 
uration space). The two stretches start and end at al- 
together four "ports" , and are connected to each other 
through two external "links" 13] . Changing the way the 
ports are connected inside the encounter, as illustrated 
by the dashed line in Fig. ^1, we obtain two disjoint pe- 
riodic orbits, each containing one link and one encounter 
stretch. The existence proof of the latter two orbits as 
exponentially close to the links of the original orbit re- 
quires chaos and was in essence given in Ref. These 
two orbits can either both be included in B such that D 
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FIG. 1: "Diagrams" of simplest (pseudo-)orbits in configura- 
tion space. Arrows indicate sense of traversal. 



remains empty, or vice versa, or one is included in B and 
the other one in D. At any rate the cumulative actions 
of (A, C) and (B, D) almost coincide. 

A second example, a so-called Sieber/Richter pair 
is shown in Fig.^>: Here two orbits differ in an encounter 
wherein two stretches are approximately mutually time- 
reversed. One of these orbits is included in A leaving C 
empty or vice versa, whereas the other is included in B or 
D. Such orbit pairs exist only in time-reversal invariant 
systems since we assume that the sense of motion on an 
encounter stretch or a link can be reverted in time. 

More complicated quadruplets involve any number of 
orbits, and (^4, C) and (B, D) can differ in any number of 
encounters where arbitrarily many stretches come close 
in phase space (modulo time reversal for time-reversal 
invariant systems); encounters may be self-encounters 
within periodic-orbit components in a pseudo-orbit, mu- 
tual encounters of different periodic orbits within one 
pseudo-orbit, or even from different pseudo-orbits. 

As in the diagonal approximation some of the compo- 
nent orbits of (A, C) may just be repeated in (B, D). To 
evaluate the generating function, we must sum over all 
quadruplets of this type. We can split the sum into one 
over the "diagonal" parts of these quadruplets and one 
over the orbits differing in encounters. The first subsum 
yields Zdiag as in JHJ such that Z = Zdi ag (l + Z g) with 

Z oS = (F A F* B F c F* D (-l) n ° +n ° (10) 

diff. in one. 

To find Z s we must classify and count quadruplets 
with encounters. Their topological "structure" must be 
dealt with first. To that end we number all encounter 
stretches (stretches, for short) of A and C. Starting with 
an arbitrary stretch in the orbit with the smallest num- 
ber of stretches, we number the stretches in that orbit 
in their order of traversal; we continue with the orbit 
with the second-smallest number of stretches, etc. Each 
structure now corresponds to one way of (i) grouping 
these numbered stretches into encounters, (ii) choosing 
their mutual orientation if the system is time-reversal in- 
variant, (iii) changing connections inside the encounters, 
and (iv) dividing the original orbits among A and C and 
the reconnected orbits among B and D. Quadruplets of 
pseudo-orbits can typically be described by several equiv- 
alent structures: If we denote by w M the number of com- 
ponent orbits of j4, C with /j, encounter stretches, there 



are 11^/^^' wavs of choosing first stretches, and or- 
dering the LUfi orbits with fixed ji. When summing over 
structures, we must divide out that overcounting factor. 

Next, pseudo-orbit quadruplets are characterized by 
phase-space separations between the encounter stretches. 
To measure separations for an /-encounter, we introduce 
a Poincare surface of section orthogonal to the original or- 
bit in an arbitrary reference point in one of the stretches. 
The other stretches pierce through the same section in 
(I — 1) further points; their phase-space separation from 
the first piercing can be decomposed into components 
Ui and Si along the unstable and stable manifolds. As 
shown in Rcf.jJ], the encounter contributes to the ac- 
tion difference as AS — SjUj and has a duration 

tone = t In rr-T3 l. i where A is the Lyapunov 



A maxj |sj|xmaxfc \uk\ 

exponent and c a constant whose value is unimportant. 

The sum over A, B, C, D in i|10[l can be written as a 
sum over structures and an integral over s, u and the 
link durations t. The measure to be used 0, 0] ob- 
tains a factor l/f2' -1 t onc from each /-encounter, with 
the volume of the energy shell. The factor arises 
from ergodicity and gives the uniform probability den- 
sity for the / — 1 later piercings to have given stable and 
unstable differences from the reference one; the factor 
1 /t e nc compensates an overcounting due to the fact that 
the Poincare section may be placed anywhere inside the 
encounter. 

We next split the phase-space integral into factors 
representing the links and the encounters. To do so, 
we write the time Ta as a sum of durations of all 
links and encounter stretches which belong to A before 
reconnection; Tg, lb, and Tr, are decomposed similarly. 
We then obtain an integral J Q °° dt c^^a o, c~ £ b or d)/ Th 
for each link (belonging to A or C before recon- 
nection and to B or D afterwards) and an integral 



fdf-hidf-hi 



for each encounter (with Ia-,Ib-,^c^d the numbers of 
stretches of the encounter belonging to A, B,C, D, and 
I A + lc = Ib + Id-) Evaluating these integrals as in 0] 
we obtain a factor i(e\ or c ~ 6 b or d) 1 f° r eacn un k a 
factor i(l A £~A ~ Ib£~b + Ic^c ~ ^£> £ l>) f° r eacn encounter, 
while Th cancels out. In this way, Z g becomes the sum 
over structures 



-"off ' 



E «"- B+ " £> Ifcnc ^A<ii-lBe-+lce++l D € D ) 
(-i)"c+" D n M ^^!n link3 (-i)(4„ r c- e io, 

struct 



(11) 

where for (3=1 the factor k" b+ " b accounts for the two 
different senses of motion on the "reconnected" orbits. 
Summarily referring to the linear combinations of the 
^A bcd m t ne link an( l encounter factors as e, we infer 
that Z s is a power series in 1/e. The term (l/e) m is 
provided by all structures with m = L — V, with V the 
number of encounters and L the number of stretches in 
a structure; note L > V. This remark allows to draw all 
"diagrams" contributing to each of the first few orders of 
the expansion and to evaluate their contributions. 
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For instance, the order m = 1 is determined by the two 
diagrams in Fig. ^ whereas for m — 2 we need quadru- 
plets with two 2-encounters or one 3-encounter. In the 
unitary case, all these quadruplets either yield vanish- 
ing contributions to Z (after summing over all possible 
assignments of orbits to A, B, C, D) or mutually cancel. 
Reassuringly, the low-order analysis of the general ex- 
pression above complies with the fact that for (3 = 2 the 
diagonal approximation exhausts the RMT result. 

In the orthogonal case, off-diagonal contributions re- 
main and the non-oscillatory and the oscillatory parts 
of the correlator C(e) are obtained according to 
Not surprisingly, the non-oscillatory terms are deter- 
mined only by pairs of orbits (such as Fig. \\}p) known 
from the previous work on the small-time form factor; 
in the present language, either A or C, and either B 
or D are empty. All genuine pseudo-orbits end up con- 
tributing nothing. The first oscillatory term, oc e 2lc /e 4 , 
does involve non-trivial pseudo-orbit quadruplets. It 
can be attributed to quadruplets of orbits involving two 
Sieber/Richter pairs; further (mutually canceling) dia- 
grams are archived in an appendix. Proceeding to all or- 
ders we get the full asymptotic expansions (JTJ a la RMT. 

The form factor is finally obtained from the complex 
correlator as KM = ^ f°° + ? a e~ a " C(e)de, a > 0. 
Substitution of our expansion Q for C(e) and term-by- 
tcrm integration produces, for both symmetry classes, 
two analytic functions different from zero for Rer > 
and Rer > 1, respectively, which sum up to the RMT 
form factor. The Laplace transform of the result repro- 
duces the exact complex correlator of RMT as C(e) = 
2 J °° e l2T€ K(t)oIt including the power law of its real 
part for e — > 0, i.e., level repulsion. For conditions un- 
der which such a back-and-forth Laplace transformation 
(Borel summation) uniquely restores a function from its 
asymptotic series see fl5j . 



We conclude with the following remarks. Implicit 
to our present analysis is a specific order of two lim- 
its. These are the semiclassical limit which brings in 
the periodic-orbit sum a la Gutzwiller in QJ, alluded 
to as liniA^o below, and the vanishing of the imagi- 
nary part lmE + = "/A/2t: of the complex energies and 
of Ime + = 7, i. e. lim 7 ^o- We need the condition 
7 A/ A = 7 > 1 to make the periodic-orbit contributions 
to our asymptotic expansions well defined. It is worth 
noting that this condition effectively limits the orbit pe- 
riods as T < Th', see Eq. (IKJf) where the final exponential 
includes a damping q-^( t a+Tb+t c +t d )/t h _ j n our lixxiit 
sequence (first A — > with 7 > I fixed, then 7 — > 0) the 
two representations, x and |j, become inequivalent, and 
resolve different parts of the two-point function (oscil- 
latory and non-oscillatory); separate asymptotic expan- 
sions for both are required to recover the full informa- 
tion. This complementarity phenomenon is reflected in 
the structure of the field theoretical approach to spectral 
statistics as well: the functional integral representation of 
Z in RMT is controlled by two saddle points . In the 
limit 7< 1, both saddles equally contribute and give C 
in full in either representation, || and x; one saddle pro- 
vides the non-oscillatory part of C, the other the oscil- 
latory part. However, for 7 > I, the oscillatory part gets 
exponentially suppressed as e -7 in the ||-representation, 
while the x -representation has the non-oscillatory part 
exponentially damped. The perturbative e~ ^expansions 
of the surviving parts, on the other hand, quantitatively 
agree with the results of our semiclassical analysis, to all 
orders in 1/e. We shall further discuss analogies semi- 
classics/field theory in a forthcoming paper. 
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Appendix: Leading contributions 

To illustrate our results, we evaluate the contributions 
of the simplest quadruplets of pseudo-orbits. We start 
with the "diagram" in Fig. la, where an orbit involving 
an encounter of two almost parallel stretches decomposes 
into two disjoint partner orbits. First, assume that the 
initial orbit is included in the pseudo-orbit A. Each of 
the two partner orbits can be included either in B or in 
D. This leads to four different possibilities, whose contri- 
butions to Z g can be evaluated using the diagrammatic 
rules introduced above, 



2ie\ — 2ie B 



- 2 



2ie+ 



H(4-%)) 2 H(4-^))H(4-^)) 



(12) 



+ 



2iet - 2ie 



D 



H(4-e B )) 2 



= 0. 



In a similar way, all possibilities with the initial or- 
bits included in B,C or D also sum to zero. The orbit 
quadruplets corresponding to the "diagram" in Fig. la 
thus yields no net contribution. 

Next, we consider Sieber/Richter pairs as in Fig. lb. 
These pairs require time-reversal invariance and there- 
fore exist only in the orthogonal case. According to our 
diagrammatic rules, Sieber/Richter pairs where, say, the 
original orbit is included in A and its partner is included 
in B yield a contribution + 2 — to Z Q g . All other pos- 

sibilities to assign the two orbits to A,B,C,D can be 
treated in the same way. Summation thus yields 2F, 
where 



F 



i(4 



(13) 

Upon taking derivatives and identifying the energies 
in a columnwise way, Sieber/Richter pairs yield a non- 
oscillatory contribution jrk^ to the complex correlator, 
whereas the corresponding oscillatory term (obtained af- 
ter crosswise identification of energies) vanishes. 



The simplest quadruplets beyond Fig. la and b are 
those with m = L — V = 2, i.e., either with two 
2— encounters or one 3— encounter. These quadruplets 
are shown in Fig. 2, grouped into columns labeled by 
1 : 1, 1 : 2, 1 : 3, and 2 : 2. Here the first number counts 
the orbits before reconnection (black lines) whereas the 
second number counts the orbits after reconnection (col- 
ored or gray lines). To evaluate the resulting contribu- 
tions, we have to sum over all ways of distributing these 
orbits among A, B, C, and D (placing the former orbits 
in A or C and the latter in B or D or vice versa) and 
apply the above diagrammatic rules. 

We have done so, and obtained the following result: All 
quadruplets of pseudo-orbits that do not require time- 
reversal invariance (in the shaded rows) give mutually 
canceling contributions to Z a s, and thus no off-diagonal 
contributions to the correlator. 

For time-reversal invariant systems Z g is non-zero. 
When taking derivatives and identifying energies in a 
columnwise way, we see that non-oscillatory terms in 
C(e) arise only from quadruplets associated to 1 : 1, sum- 
ming Up tO 2(^jT- 

In contrast, non-vanishing oscillatory contributions 
(obtained through the crosswise representation) arise 
only from the quadruplets of type 2:2. The first and 
third of these diagrams in Fig. 2 mutually cancel. The 
leading off-diagonal contribution to the oscillatory part 
can thus be attributed solely to the second diagram, com- 
posed of two Sieber/Richter pairs . For these quadruplets 
the sum over possible assignments to A, B, C, and D 
can be performed independently for the two pairs. We 
are therefore left with the square of the contribution of 
a single Sieber/Richter pair, which must be divided by 2 
since the divisor I1a m TO ^ s now ec l ua l to 2. We 



) thus obtain a term 2F 2 in Z g and a contribution 
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to C(e), coinciding with the leading oscillatory term in 
the RMT prediction (JTJ. 





FIG. 2: "Diagrams" responsible for the asymptotic expansion of the complex correlator up to order O(-j). The shaded 
rows contain pairs which exist even in the absence of time-reversal invariance; the rest exists only in the orthogonal case. 
The quadruplets in the first column are marked by 1 : 1 since they involve one orbit before and after reconnection. In 
others columns, the number of orbits before and after reconnection is 1 : 2, 1 : 3 and 2 : 2, respectively. Pseudo-orbits 
containing a 3-encounter can have several partners depending on the reconnection. The contributions to Z a s as defined in 
llll l are shown explicitly for pairs with L — 2 links and V = 1 encounter. For m = L — V = 2, the resulting contributions 
to Z f{ are rather involved and will be presented in a forthcoming publication. 



